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ABSTRACT 

We present critical, long-wavelength observations of r\ Carinae in the submillimetre using 
SCUBA on the JCMT at 850 and 450 /im to confirm the presence of a large mass of warm dust 
around the central star. We fit a two-component blackbody to the IR-submm spectral energy 
distribution (SED) and estimate between 0.3 — 0.7 M© of dust exists in the nebula depending 
on the dust absorption properties and the extent of contamination from free-free emission at 
the SCUBA wavelengths. These results provide further evidence that 77 Carinae's circumstellar 
nebula contains > 10 M Q of gas, although this may have been ejected on a longer timescale 
than previously thought. 

Key words: circumstellar matter-stars: individual: rj Carinae-submillimetre 
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1 INTRODUCTION 

rj Carinae is one of the most luminous objects in our Galaxy 
(10 6 Lq) and is well known for dramatic outbursts in which ma- 
terial is ejected outwards from the star. Its most famous mass ejec- 
tion is thought to have occurred in the early 19th Century, where a 
few solar masses of gas was expelled in a few decades. This feature 
is clearly seen in observations as a bipolar nebula, known as the 
Homunculus (Gaviola 1946). The origin of the shape of the Ho- 
munculus is unknown, but could occur due to ejection of mass in 
a non-uniform interstellar medium i.e. if the ambient medium was 
denser in the equatorial plane, the nebula would preferentially ex- 
pand along the poles. Surrounding the Homunculus is a larger, less 
dense nebula, known as the 'outer ejecta' (e.g. Weis 2005) which 
contains numerous filaments and dense condensates not dissimilar 
to those seen around the young supernova remnant Cassiopeia A 
(e.g. Smith & Morse 2004, Fesen et al 2001). The central star is 
hidden from view by this structure but is thought to be a luminous 
blue variable (LBV) star with mass 100 M (Pittard 1999; Smith 
et al. 2003). LBV stars are well known for large luminosities, vio- 
lent instabilities and periods of large mass loss. The presence of a 
binary companion (e.g. Pittard 2003) would provide an explanation 
for the catastrophic mass ejections and the strong variabilities seen 
in the radio and X-ray (Cox et al. 1995; Duncan & White 2003). 
The cycle of periodicity suggests that any exisiting binary would 
have a highly eccentric orbit, leading to strong wind-wind colli- 
sions such as those seen in WR binaries (Marchenko et al. 2002; 
Pittard 1999). The evolutionary phase of an LBV star can last more 
than 10 4 yrs and hence in this phase alone, a single star could lose 
roughly 1 Mq of gas. Episodic eruptions will increase the mass loss 
further. It is therefore not unreasonable to assume that the dust mass 



loss rate could be > 10" 6 M 
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For the lifetime of the aver- 



age LBV phase, we could expect LBV stars to inject ^ 0.01 Mq 
of dust into the interstellar medium and could be responsible for 
the dust emission seen around the young supernova remnants Cas 
A (Dunne et al. 2003, Wilson & Bartla 2005), Kepler (Morgan et 
al. 2003), SN1998J (Pozzo et al. 2004), the Crab (Green, Tuffs & 
Pospescu 2004), SN2002 (Barlow et al. 2005) and more recently 
SN2003gd (Sugerman et al. 2006). 

The structure of rj Car is revealed in the many multiwave- 
length observations, yet there are issues surrounding the origin of 
the Homunculus. IR observations at 12 - 17 /im revealed a small 
structure within the Homunculus which appeared to be concen- 
trated around the 'waist' of the nebula (Morris et al. 1999; here- 
after M99). They found a small bright core a few arcsecs across 
which they interpreted as a disk/torus of dust with two components 
at 200 K and 1 10 K with ~ 0.15 M Q of dust. They suggested (with 
a normal gas-to-dust ratio of 100) that the gas within this small 
torus could therefore be ~ 15 Mq . This argument is flawed since 
estimating the 1 10 K dust mass relied on observations at the longer 
infrared wavelengths where resolution is limited to greater than 1 
arcminute, so the precise location of the 1 10K dust component can- 
not be determined from the mid-IR images alone. 

Davidson & Smith (2000) also argued against the IR emis- 
sion in a small torus, proposing instead that this was too small to 
account for the large luminosities radiated at 100 fim. The visual 
image provides further evidence for this argument since there is no 
evidence of an opaque dusty disk on such a small scale. Hony et 
al. (2001) then used high resolution images at 8 - 20 /im to suggest 
that the observed IR dust is not warm dust at 1 10 K, but actually hot 
dust at 200 K with a strange overlapping ring structure in the mid- 
dle of the nebula. They argue that these rings are similar to those 
seen in SN1987A and are at different polar axes to each other. If 
this scenario is correct, the total dust mass in 77 Car decreases to 
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Figure 1. r] Car at 450 fim with 850 fim contours overlaid. White contours 
start at 5-cr (160 mJy/beam) incrementing by 2-cr with black contours start- 
ing at 10-O-. Both images have been smoothed to the same resolution and 
are in J2000 coordinates. 



and sky dips to give relatively stable optical depths of T450 ~ 0.8 
and T850 ~ 0.2 (see Section ITHl . The remsky option was used 
with the outer ring of bolometers flagged as sky bolometers (as 77 
Car is almost a central point source in the maps) and the option for 
addback was turned on. 

The 450 fim SCUBA image of rf Car is shown in Fig. \\\ 
with 850 crni contours overlaid. White contours start at 5-cr (160 
mJy/beam) incrementing by 2-cr with black contours starting at 10- 
a. The central part is not resolved by the SCUBA field-of-view but 
we see an interesting feature which appears to be extended in the 
x direction, along the midplane of the image which is particularly 
clear at 850 fim. It is difficult to determine the significance of this 
since the source is at such a high airmass. It is well known that 
larger chop throws can cause 'smearing' along the direction of the 
throw, in this case the x-axis, and the image of 77 Car is similar 
to example images of this problem 1 . Extended emission along the 
midplane is also seen in the calibrator CRL618 at 850 fim (Figs [2] 
(b) & (d)), suggesting it may be a combination of the chopping 
method and the beam pattern at high airmass. The extension is not 
as clear in the 450 mum image of the calibrator. Radial profiles of 
both the beam (taken from the calibrator) and the source show that 
the extension cannot be separated from the beam pattern. 



~ 10~ 4 Mq due to the strong dependence of dust mass with tem- 
perature; the gas mass ejected in the Homonculus would then be 
< 10 -2 Mq and the evolution of 77 Car would not be as violent 
as currently thought. This claim has been disputed in Smith et al. 
(2000; 2003) who propose that the strange smooth loop structure 
seen in the Hony et al. observations are simply due to a disrupted 
torus. They point to compact knots and hotspots in the IR with fila- 
ments and arcs which could be due to a toroidal distribution of dust 
which has been disrupted by subsequent mass loss. They also use 
temperature maps from ISO data to show that cooler 'hidden' dust 
may exist in the outer polar lobes of the Homunculus and not in a 
compact torus. 

M99, Davidson & Smith (2000) & Smith et al. (2002) have 
all argued that there is missing mass which has never been ac- 
counted for by the observations, and the cooler dust mass which 
has been found has not been resolved, so its precise location is 
unknown. If there is missing ejected mass, this will have impor- 
tant implications for the energy budget and estimated parameters 
of LBV stars as well as its future evolution. To determine if colder 
dust could exist further out from the centre of the star, wavelength 
measurements are needed beyond the mid-IR range (Davidson & 
Smith 2000) which will dominate the dust mass. We therefore use 
archived SCUBA observations of 77 Car to investigate the distribu- 
tion of the dust and the ejected mass in one of the most studied 
objects in the nearby Universe. 



2 SUBMILLIMETRE OBSERVATIONS 

The sub-mm observations of 77 Carinae were originally observed 
by SCUBA on the JCMT in January 1998 in Grade II weather, sim- 
ulataneously observing at 450 and 850 fim (with nominal beam 
size of 8 and 14 arcsec respectively). The data was taken from the 
JCMT archive (originally observed by Henry Matthews, program 
ID M97bc30) and consisted of a single Jiggle map observation 
with a 120 arcsec chop throw in Alt-Az and position angle 90°. 
This was reduced using the standard data reduction pipeline SURF 
(Sandell et al. 2001) with sky opacities obtained from CSO-FITS 



3 CALIBRATION & ERRORS 

The high airmass of 77 Carinae (~ 5.5) at the time of observation 
produces potential difficulties with calibration which we will ad- 
dress carefully in this section. The calibrator observed on the same 
night was at a lower airmass of 2.5. Calibration of the data requires 
two steps, firstly an extinction correction for the absorption of in- 
coming radiation by the atmosphere, and secondly a comparison 
with a calibrator of known flux to convert instrumental units (V) 
into Jy. 

3.1 Extinction correction: 

For the time the data was taken, the JCMT has two methods avail- 
able to measure the value of the zenith sky opacity (r) at the wave- 
lengths of interest; skydips and extrapolation from the CSO 1.1 mm 
t. Skydips are taken at the SCUBA wavelengths (450 and 850 fim) 
but are only made infrequently (every two hours or so). The CSO 
radiometer measures the r at 1.1 mm every 10 minutes, and so sam- 
ples changes in the atmosphere on a much more useful timescale. 
The JCMT have created relationships between the CSO r values 
and those at 450 and 850 fim using an extensive data-base of sky- 
dips. They also fit the CSO r data for each night with a polynomial 
which allows interpolation between measurements and reduces the 
effect of any 'spikes' in the CSO measurements. These polynomial 
fits are the recommended way to extinction-correct SCUBA maps. 

On the night the data were taken, the sky was extremely stable. 
The plot of the CSO r values, the polynomial fit and the residuals 
is shown in Fig. [3] The part of the night considered to be 'stable' 
after sunset and before sunrise (when the 77 Car data were taken) is 
shown between the dashed lines. 

The value of r which is used to correct the flux of the object 
for atmospheric absorption, is uncertain - and this uncertainty, ov, 
leads to an uncertainty on the flux measured. Due to the high air- 
mass of 77 Car, any error in the r value is amplified greatly when 

1 see http://www.jach.hawaii.edu/JCMT/continuum/ 
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Figure 2. (a) Contour plot of rj Carinae at 850 jim displaying 1 percent of the peak voltage up to 10 percent, (b) Contour levels for the calibrator CRL618 
showing the same levels. Extension along the midplane is seen in both calibrator and source, (c) rj Carinae at 450 /tin. Contours starting at 15 percent of the 
peak voltage up to 95 percent, (d) CRL618 at 450 fim. Contour levels 5 percent of the peak voltage up to 50 percent. 



applied. From Dunne & Eales (2001), we have the fractional error 
in flux due to the uncertainty in r as: 



-Aa T 



(l) 



where we have taken ay to be the r.m.s. of the residuals in 
Fig. |3| and A is the airmass. The value of ay is 0.00253 and this 
must be translated to an r.m.s. on the r values at 450 and 850 /im. 
For this we use the pre-upgrade revised narrow-band filter rela- 
tionships between CSO tau and 450/850 given on the JCMT 
webpages 2 , and using propagation of errors we find that ay (850) = 
0.005 and a T (450) = 0.012. 

Using Eqn.Qwe calculate the fractional error on the flux due 
to uncertainties in the extinction corrections as a ex t (850) = 3 per 
cent and a ext (450) = 6 per cent. If we instead take the extreme 
values of the residuals rather than the r.m.s. we get increased values 
of 6 (850) and 15 per cent (450), however, these are not statistically 
correct, and we include them as merely illustrative of a worst case. 
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3.2 Changes in the gain: 

The second part of the calibration procedure is to compare the sig- 
nal measured in an aperture on the object with that in the same 
aperture on a calibrator of known flux. The calibrator used was 
CRL618, chosen because it was observed shortly before rj Cari- 
nae, and therefore sky conditions and the shape of the dish would 
be closest to those at the time of the rj Car observation. A change 
in the dish shape (due to thermal relaxation after sunset), or a tele- 
scope which is poorly focussed will result in changes in beam shape 
and hence the gain - defined as the ratio of calibrator flux (Jy) / cal- 
ibrator flux (V). Thus, if the calibrator is observed under as similar 
conditions as possible to the objects then this effect is minimised. 

A check of all of the calibration maps during the night showed 
that the beam shape was very poor in the earlier part of the evening 
(after sunset), despite many attempts at focussing. However, for the 
calibrator map taken prior to r\ Car, alignments were also made on 
the secondary mirror in the x-y plane which resulted in a much im- 
proved beam shape. Thus we believe that the previous poor quality 
images were due to a poorly aligned mirror and there is no reason 
to believe that this will affect the source data which was taken just 
after the re-alignment. 

There could be some concern that observing at such low ele- 
vations may also affect the shape of the dish and hence the gains, 

2 http://www.jach.hawaii.edu/JCMT/continuum/calibration/atmos/tau.html 



Figure 3. The CSO tau, polynomial fit and 850 jim skydips for UT 
24/01/1998. The residuals are shown below. The times used to derive the 
uncertainty in r are shown by dashed lines, rj Car was observed at ~ 0.5 on 
the horizontal axis. The larger scatter in r values at UT fractions of 0.3 and 
0.8 correspond to sunset and sunrise, where the sky becomes more unstable. 
As both our calibrator and object were observed shortly after midnight, we 
have only taken the stable portion of the night to estimate ay . 

however, the staff at the JCMT have found no evidence for a link 
between gain and elevation (Iain Coulson, private communication). 
Furthermore, the beam shapes look similar between the calibra- 
tor and rj Car (Fig.|2j, again suggesting that there is not a serious 
change in images observed at the lower elevations. Any deteriora- 
tion in dish shape tends to lead to higher gains (more Jy per V) and 
therefore any elevation effect on the dish would be leading us to 
underestimate the flux of rj Car rather than overestimating it (as the 
calibrator was observed at a lower airmass). 

Since the earlier calibration maps had been affected by the 
mis-alignment of the secondary mirror, and therefore do not reflect 
usual nightly variations in the gain, we will use the information 
from Tables 1 & 2 in Dunne & Eales (2001) to estimate the uncer- 
tainty on the gain. Here the variation of gain values in various sized 
apertures was calculated using data taken over the same period as 
these observations (July 1997 - July 1998). This gives the uncer- 
tainty in the gain in a 45 arcsec aperture as a ga in = 6 per cent at 
85 Ojim and a ga in ~ 9 per cent at 450 jim. 

There will also be an absolute calibration uncertainty due to 
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errors in the flux of the calibrator and the model used to generate 
the planet fluxes (on which the flux of the secondary calibrators is 
based). The errors on the flux of CRL618 are quoted as 4.7 ± 0.12 
Jy at 850 /xm and 11.8 ± 0.6 Jy at 450 /jm. Adding in quadrature 
an additional 5% error due to the uncertainty in the model, the total 
error on the calibrator flux is estimated as a a t s = 5.6 per cent at 
850 /im and o a bs = 7.0 per cent at 450 /im. 

The total calibration error is then the quadratic sum of the 
three terms for extinction, gain and absolute as follows: 

a cai = (J 2 ext + <J 2 galn + v 2 abs (2) 

giving (j ca ;(850) = 9 per cent and <r ca ;(450) = 13 per cent. While 
these numbers may look surprisingly low, it must be remembered 
that this was a night with exceptionally good observing conditions 
(during the El Nino period), and that the best practise has been used 
in the calibration procedure. 



3.3 The SCUBA Fluxes 

The integrated flux was measured by placing the same aperture over 
both r\ Car and CRL618, measuring the total flux in Volts on both. 
The fluxes of r\ Car at 850 and 450 /im were found to be 12.8 ± 
0.8 Jy and 52 ± 6.8 Jy respectively as measured in an aperture of 
diameter 45" . 

The radius of the aperture needed to incorporate the entire 
emission is ~ 28", which at the distance of r\ Car corresponds 
to a physical size of 0.3 pc and is just larger than the visual size 
of the Homunculus. The dust seen in these images may be probing 
the 'outer ejecta' region (e.g. Smith & Morse 2004, Weis 2004). 
However, we cannot resolve the inner 10" of the Homunculus and 
we refer the reader to the recent sub-arcsecond study of r\ Car e.g. 
Chesneau et al (2005). 



4 THE IR-RADIO SPECTRAL ENERGY DISTRIBUTION 

Due to the high variability in the X-ray, Radio, mm and IR and as 
the sub-mm observations of r\ Car were taken during a shell-event 
in 1998, we chose literature fluxes taken nearest this epoch (corre- 
sponding to a radio quiet stage). This allows for the best possible 
direct comparison of the sub-mm with any variability in the radio, 
as well as determining the evolution of the SED during this time 
(e.g. Abraham & Damineli 1999; Abraham et al. 2005). Where un- 
available, older literature fluxes were used (Cox et al. 1995). The 
1992 millimetre fluxes taken from this work are thought to be at 
a similar level to those expected in 1998 (e.g. Fig 1, White et al. 
2004). 

We fitted the SED of r\ Car from 10 - 850 /im using a mod- 
ified blackbody (Fig|4j requiring a two temperature component at 
around 200 and 110K. The dust emissivity exponent /?, is < 1.5, 
lower than the nominal value of 2 for 'normal' interstellar composi- 
tion but higher than those estimated in M99 and Smith et al (2003). 
This indicates that the grains could be amorphous, large or non 
spherical in composition (e.g. Bohren & Huffmann 1983). How- 
ever, there is likely to be contamination from the radio emission at 
the sub-mm wavelengths which is dependant on the phase of the 
radio cycle. There are two sources of radio emission from the star - 
free free ionised emission from the stellar wind and emission from 
the optically thin Homunculus region (Cox et al. 1995). Both ther- 
mal processes depend on frequency, free-free varies as u 6 (Wright 



& Barlow 1975) and the optically thin region as v~ 01 . The con- 
tamination from the radio to the sub-mm emission can be removed 
by determining the flux expected at 450 and 850 /im. The varia- 
tion of free-free emission in Jy with wavelength from an ionised 
stellar wind is given in Lamers & Cassinelli (1999). X-ray mod- 
eling of colliding wind binaries fit the observations of r\ Car with 
parameters: mass ~ 70 Mq, T c h ~ 20, 000K, mass loss rate, 
M ~ 9 x 10~ 5 M Q yr" 1 , terminal velocity ~ 9300 km s _1 and 
distance 2.6 kpc (Pittard 1999). The SED of the free-free emis- 
sion assuming these values and a doubly ionised wind is shown 
in Fig. [4] (dashed line). The power law is significantly below the 
sub-mm fluxes and also underpredicts the observed 1.2 and 3 mm 
fluxes. Cox et al. (1995) suggest that the power law slope from the 
free-free emission is actually v in the mm regime due to hydro- 
gen recombination changing the electron and proton densities in the 
wind (see their Fig 3). The emission then turns over at A ~ 1.2 mm 
to the normal slope of 0.6, just below the observed flux at this wave- 
length, shown on Fig|4| We use the observed fluxes rather than the 
theoretical prescription of a stellar wind to estimate the maximum 
possible contribution from the ionised stellar wind in the sub-mm 
regime. Thus we would expect a flux of ~ 12 Jy at 850 and 18 Jy at 
450 /im respectively due to the ionised stellar wind. It is clear from 
Fig |4| that almost all of the 850 /im flux could be free-free emis- 
sion, yet not all of the 450 /im emission can be explained from the 
expected power law. Indeed, it is difficult to explain the 450 and 
850 /im flux ratio with a radio emission process. However, given 
the uncertainty of how much of the sub-mm flux is due to dust 
emission or contaminated by free-free, we have two scenarios for 
the dust mass in 77 Car (1) best case: all of the estimated SCUBA 
flux is due to dust and (2) worst case: 94 and 35 per cent of the 850 
and 450 /im flux respectively is due to the ionised stellar wind. This 
effectively gives a lower limit to the range of dust masses expected 
(see Section|5j. 

As noted in Smith et al. (2003), the fit to the energy distribu- 
tion is not unique, a range of values of temperature could achieve 
a satisfactory fit to the spectrum within the errors on the fluxes. 
Indeed in their work, they reproduce an adequate fit to the Mor- 
ris et al. (1999) data with dust temperatures of 140 K and 200 K. 
However, none of the previous best-fit SEDs produce the correct 
450/850 ratio as observed here. The extra fluxes at the long wave- 
length tail provide another constraint, which forces the fit to have 
a colder temperature than the previous models at 110K. It is ex- 
tremely important to separate the dust emission from free-free at 
these wavlengths as this affects the dust mass in the nebula. If the 
wind is optically thin, for example, most of the IR emission would 
be due to dust and not contaminated by free-free. 



5 A LARGE MASS OF DUST? 

The emission seen in the SCUBA images and the radial profiles 
suggests that some of the dust may be located outside the Ho- 
munculus although due to the beam pattern we cannot determine 
if this is a real extension. The HST/WFPC2 image is shown in 
Fig.|5|with the CHANDRA X-ray observations and the submm emis- 
sion shown as contours starting at S/N > 3-a (original image from 
Smith & Morse 2004). The 450 /im emission appears to coincide 
with the X-ray emission which is thought to occur where previous 
mass loss from the star is colliding with the surrounding interstel- 
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Figure 4. The spectral energy distribution from 12 - 850 fim through to the radio regime of r) Carinae for the 1998 epoch. The 175 fim data point is also 
included although this was observed in 1977 (Harvey, Hoffman & Campbell 1978). The solid line is the two temperature component fit to the infra-red 
fluxes. The dashed line shows the SED expected from theoretical free-free emission (Lamers & Cassinelli 1999) with wind and stellar parameters taken from 
Humphreys (1988) and Pittard (1999). This falls short of the observed fluxes at the mm wavelengths from Cox et al. (1995). The dot-dashed line represents 
the emission from the optically thin, ionised Homunculus with spectral slope a = — O.f normalised to the 3 cm flux (Cox et al. 1995). 
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Figure 5. HST/WFPC2 image taken from Smith & Morse (2004) with con- 
tours of soft X-ray emission (thin black lines) and smoothed 450 (im (thick 
black lines). The submm contours start at 3-cr and extend just outside the 
Homunuculus region shown by the well known HST bi-polar nebula and 
may be associated with the outer ejecta. Indeed there is a good correlation 
between the location of the 450 fim emission and the hard X-Rays from the 
shocked wind. The image was aligned using the position of the central star. 
Co-ordinates are J2000. 



lar material, causing shock heating of the gas. We may therefore be 
probing dust shells from a previous mass-loss phase 3 . 



3 A rough calculation shows that the dust could not be due to swept up 
from interstellar gas since the swept up dust mass at this distance would be 



To determine the dust mass in r\ Car, we use 
S V D 2 



k v B v {v,T) 



(3) 



less than < 10~ J M 



0- 



where k„ is the dust mass absorption coefficient (~ 0.27 kg -1 m 2 
for normal interstellar dust at A45o,jm) and D is the distance to the 
star. The dust mass from the best fit SED to the free-free subtracted 
fluxes (scenario 2) is < 0.54 ±0.1 Mq depending on the absorp- 
tion properties of the dust. With an typical gas-to-dust ratio of 100, 
this suggests that around 50 Mq would have been expelled in the 
wind. Such a large number may indicate that the gas-to-dust ratio 
in the atmospheres of massive stars could be higher than predicted 
(e.g. AG Car - Voors et al. 2000). Alternatively, the dust could be 
more efficient at emitting than normal ISM dust. Previous mod- 
elling of the SED of 77 Car has shown that the dust does not follow 
an opacity law of k a A _/3 (Robinson et al 1987). If we use the 
absorption coefficient estimated from Dunne et al. (2003) to repre- 
sent pristine, newly formed dust in stellar winds, we obtain 0.3 Mq 
which would require ~ 30 Mq of gas to be expelled. Using the best 
fit fluxes without the free-free contamination (scenario 1), the dust 
mass would be 0.7 ± 0. 1 Mq . Errors quoted are given using 68 per 
cent confidence intervals from parameters obtained using the boot- 
strap technique (whereby 1000 sets of artificial fluxes were created 
from the original data with each individual set then fitted with our 
two- temperature model). 

This is in good agreement with the original hypothesis of M99 
and Smith et al. (2003). The dust mass calculated here appears to 
be spread out over a larger area than the torus and quite possibly 
outside the polar lobes of the Homunculus. Combining this result 
with that in M99 and Smith et al. (2003), we obtain the dust masses 
for each component as listed in TableQ Since we cannot resolve the 
inner Homunculus, the dust mass here in the outer ejecta is likely 
to include the mass postulated to be in the lobes or the torus. 

If we assume that the dust around rj Car is similar to dust 
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Region 


Size (") 


Temperature 


Dust Mass (M ) 


Core 


0.3 


400 


3 X 10~ 4 * 


Lobes, torus 


5 


200* 


0.02* 


Lobes 


20 


140* 


0.13* 


Lobes & Outer ejecta? 


28 


114 


0.54 ±0.1** 



Table 1. Components fitted to the IR-submm SED at various temperatures. 
'Average results published in M99 and Smith et al. (2003) are included for 
completion. They used a dust emissivity of 1 - 1.2 in their work. ** Dust 
mass obtained using the best-fit parameters to the SED with the free-free 
subtracted SCUBA fluxes, scenario (2). Errors are 68 per cent confidence 
intervals estimated using the bootstrap technique. 



in the ISM, then the Homunculus and the surrounding area con- 
tains roughly 0.5, M which has been ejected during the the last 
10 4 yrs. The averaged dust mass loss rate over the entire phase is 
then ~5x 1CP 5 M© yr _1 . This is an order of magnitude greater 
than the dust mass loss rate estimated from observations of Wolf- 
Rayet binaries (Marchencko et al. 2002) and more than suggested 
by previous IR observations of r\ Car. If the submm dust is located 
in the polar lobes of the nebula then it must have been ejected dur- 
ing the same singular event which caused the lobes - the great erup- 
tion. This eruption lasted only a few decades which gives a large 
time-averaged gas mass loss rate of > 0.5 Mq yr _1 (for more ef- 
ficient dust emitters). Such a high mass loss rate was first noted in 
Smith et al. (2003) and suggests that the pulsational mass loss in r? 
Car is far more extraordinary than imagined. 

The sub-mm observations suggest that an upper limit of 
0.7 Mq of dust could exist in the extended dust shells, dust which 
must have been formed in the stellar wind with effective temper- 
atures T c ff ^> 2000 K. This indicates that contrary to the study 
in Morgan & Edmunds (2003), significant amounts of dust could 
form in the stellar winds of massive stars even when the conditions 
in the atmospheric envelope are not reached by the available hy- 
drodynamical models. This previously undetected mass of material 
in the outer ejecta has important consequences for our understand- 
ing of stellar evolution and the formation of dust in massive stars. 
Massive stars may therefore provide a sustainable dust yield to the 
ISM, in which case their contribution in chemical evolution models 
could be severly underestimated (e.g. Whittet 2003). 

We have used critical long wavelength observations to show 
that a large mass of dust does exist around r\ Carinae and sug- 
gest that the dust is distributed over a larger region than previously 
thought. This has important implications for stellar evolution and 
may provide a natural explanation for the origin of the dust seen in 
young supernova remnants. However r/ Car is extremely massive 
and highly unstable, so the huge dust mass loss rate is not unex- 
pected. It is very unlikely that the lower mass stars which were 
responsible for the supernova events of Kepler and Cas A would 
have ejected enough mass to form such a large quantity of dust 
before the explosion. In order to determine a more accurate dust 
mass in r\ Car, we require simultaneous observations in the sub-mm 
and mm to separate the contamination of the free-free emission. 
We also require better resolution to determine how extended the 
sub-mm emission is. Future observations with Herschel and SMA 
will allow us to overcome this problem. These observations should 
also allow us to determine any variability of the star in the submm. 
Quantifying the dust injection from stars of different masses in dif- 
ferent stages of their evolution will have important consequences 
for estimating the dust budget in both local Galaxies and those at 



high redshifts as well as providing fundamental information about 
stellar evolution. 
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